rnormsTest <- function(){
	N<-1000	#入力次元
	XRandSd <- 1/N	#標準偏差
	XRandM <- 0	#平均
	r<-matrix(rnorm(N,XRandM,XRandSd),N)
	hist(r)
	par(new=T);curve(dnorm(x,XRandM,XRandSd),-3/N,3/N)
	cat(sprintf("ガウス乱数ベクトルxの2乗平均<x^2>の誤差率=1 … %f
	",sum(r*r)*N))
	cat(sprintf("ガウス乱数ベクトルxの分散の誤差率var(x)*N=1 … %f
	",var(r)*N*N))
	cat(sprintf("ガウス乱数ベクトルxのノルム=1 … %f
	",sqrt(sum(r*r)*N)))
}
rnormsTest()
